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Abstract. We have observed 60 edge-on galaxies in the 
NIR in order to study the stellar distribution in galax- 
ies with box/peanut-shaped bulges. The much smaller 
amount of dust extinction at these wavelengths allows us 
to identify in almost all target galaxies with box/peanut- 
shaped bulges an additional thin, central component in 
cuts parallel to the major axis. This structure can be iden- 
tified with a bar. 

The length of this structure scaled by the length of the 
bulge correlates with the morphologically classified shape 
of the bulge. This newly established correlation is there- 
fore mainly interpreted as the projection of the bar at 
different aspect angles. Galaxies with peanut bulges have 
a bar seen nearly edge-on and the ratio of bar length to 
thickness, 14 ± 4, can be directly measured for the first 
time. In addition, the correlation of the boxiness of bulges 
with the bar strength indicates that the bar character- 
istic could partly explain differences in the bulge shape. 
Furthermore, a new size relation between the box/peanut 
structure and the central bulge is found. Our observations 
are discussed in comparison to a N-body simulation for 
barred galaxies (Pfenniger & Friedli 1991). We conclude 
that the inner region of barred disk galaxies are build up 
by three distinct components: the spheroidal bulge, a thin 
bar, and a b/p structure most likely representing the thick 
part of the bar. 

Key vifords: Galaxies: evolution - Galaxies: spiral - 
Galaxies: statistics - Galaxies: structure 



1. Introduction 

Our statistical study of ~ 1350 edge-on galaxies derived 
from the RC3 ( Third Reference Catalogue of Bright Galax- 
ies, de Vaucouleurs et al. 1991 ) (Liitticke et al. 2000a, 
hereafter Paper I) has revealed that 45 % of all bulges (SO- 
Sd galaxies) are box- and peanut-shaped (b/p). This high 
frequency must be considered in a theory of bulge forma- 



tion and especially in scenarios explaining the formation 
process of b/p bulges. 

Currently discussed theories of bulge formation are 
primordial or secular evolution scenarios. While mono- 
lithic collapse (Eggen et al. 1962 (ELS)), clumpy collapse 
(Kauffmann et al. 1993, Baugh et al. 1996| ), and inside-out 
formation (van den Bosch 1994 Kepner 1999) are pro- 
posed for primordial bulge formation, merging processes 
(Wyse et al. 1997 and references therein) and dynamical 
evolution triggered by bars (e.g. Combes et al. 1990) are 
the basis of secular evolution scenarios. However, combi- 
nations of both scenarios are also possible (Combes ^000 ). 

In the last 20 years several working hypothesis about 
the origin of b/p bulges have been developed. External 
cylindrically symmetric torques can act on an initially 
non-rotating triaxial spheroidal stable system so that it 
becomes box-shaped (May et al. 1985 ). However, such 
torques are "problematical to rel ate to some astrophys- 
ical counterparts" (Combes et al. 199C). Mergers of two 
disk galaxies are proposed by Binney & Petrou ( 1985 ) 
and Rowley ( 1988 ) as a possible formation process of b/p 
bulges. However, such a scenario requires very special con- 
ditions such as the precise alignment of the spin and or- 
bital angular momenta of the two galaxies (Bureau 1998 ). 
A more common process to form b/p bulges could be the 
accretion of satellite galaxies, although an oblique impact 



angle of the sate llite is needed (Binney & Petrou 1985 
Whitmore 



* Based on observations collected at ESO/La Silla (61. A- 
0143), DSAZ/Calar Alto, and TIRGO/Gornergrat. 



Bell |1988D and the mass of the accreted satel- 
lite is constrained. If it is too large, the stellar disk would 
be disrupted (Barnes |1992| , Hernquist |1993D . The origm 
of b/p bulges by bars was first proposed by Combes & 
Sanders (1981). This theory was strongly supported ten 
years later by the studies of Combes et al. ( 19901 ), Raha 
et al. (1991), and Pfenniger & Friedli (1991) dealing with 
N-body simulation of disk galaxies. 

The connection between b/p bulges and the presence 
of a bar is difficult to prove by surface photometry, be- 
cause b/p bulges are observable only in almost edge-on 
galaxies. However, there is photometric evidence in a few 
edge-on galaxies from cuts parallel to the major axis (de 
Carvalho & da Costa 1987, Dettmar & Ferrara 1996) and 
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in two intermediately inclined galaxies (Bettoni & Gal- 
letta 1994, Quillen et al. 1997) pointing to bars. Addition- 
ally, kinematic observations of around 20 b/p bulges with 
characteristic bar signatures in velocity fields and in orbit 



analyses (e.g. Kuijken & Merrifield 1995) are reported. In 
addition, our statistics (Paper I) have shown that the over- 
all frequency of barred galaxies (~ 55 %, Paper I, Knapen 
et al. 2000) is able to explain the high fraction of b/p 
bulges. 

On the other side observational evidence for accreted 
material in b/p bulges is found by Bureau & Freeman 
( |1999D for a few b/p bulges. Indirect evidences are sub- 
structures and asymmetries in some special b/p bulges 
(Dett mar [1989| , Dettmar & Barteldrees |1990a| Liitticke 
et al. pOOObj, he reafter Paper III). Combes et al. ( [1990D 
and Shaw ( 1993| ) conclude from their investigation of b/p 
bulges , listed by J arvis ( 1986| ), de Souza & dos Anjos 
( 1987 ), and Shaw ( 1993 ), that they show less evidence 
for significant merging and are likely not the result of re- 
cent accretion events. However, b/p bulges can indirectly 
originate from interaction and finally accretion by trig- 
gering the formatio n of b ars (Noguchi 1987, Gerin et al. 
19901 , Walker et al. 199(]| ), whic h the n actually b uild u p 



the b/p structures (Fisher et al. |l994 Mihos et al. [19951 ). 
Therefore the search for bars is one of the best ways to 
discriminate between the proposed evolutionary scenarios 
of b/p bulges. 

The essential advantage of investigations in the NIR is 
the small extinction by dust at these wavelengths (Rieke 
& Lebowsky 1985 , Knapen et al. 1991). Therefore the in- 
fluence of the dust lane in the plane of edge-on galaxies 
is largely reduced by NIR observations, especially in the 
-fC-band. In this way the structure of a galaxy near the 
plane in radial as well as in vertical direction can be anal- 
ysed and thereby bars detected (Dettmar & Ferrara 1996| ) . 
In addition it is even possible to determine parameters of 
bars, e.g. their thickness and projected length. Therefore 
it is possible to investigate the connection between bars 
and shapes of bulges. 

2. Observation and data reduction 

Observations were carried out at the 1.5m TIRGO/Gor- 
nergrat in 1994/95 (run identification: TI) with NICM0S3 
ARNICA (0.92'7pixel), the 1.2m telescope at Calar Alto 
in June 1996 (CA4) and June 1997 (CA2) with MAGIC 
(Rockwell NICM0S3 array) using the high resolution 
mode (1.2"/pixel), the 2.2m telescope at Calar Alto during 
two runs in January 1997 (CA3) and March 1998 (CAl) 
also equipped with MAGIC (hr mode: 0.64"/pixel), and 
at ESO/La Silla with the 2.2m telescope in May 1998 (El) 
with IRAC-2b (NICM0S3 array) using the LC objective 
(0.49"/pixel). To the resulting sample of 56 galaxies four 
galaxies kindly provided by Emsellem (E2: ESO/La Silla 
2.2m telescope, April 1998, IRAC-2b, LC objective) are 
added. The sample was selected from the list of edge-on 



galaxies with bulge classification (Paper I) according to 
the allocated observing time, matching the field size, and 
covering all bulge types. Filters and typical on-source in- 
tegration times for each observing run are listed in Table |^ 
and the final sample is given in Table ^j. At ESO /La Silla 
and TIRGO the K'- resp. if-band filters were used. Due 
to large thermal emission in K for MAGIC we observed 
at Calar Alto mainly in H. Using different filters we can 
show that the bulge shape in the NIR is independent from 
the observed wavelengths (Paper I). 

Table 1. Observational data of the NIR 



(1) 


(2) 


(3) 


(4) 


(5) 


obs. run 


f 


int. time 


seeing 


num. of gal. 


TI 


K 


40 


1.8" 


10 


CAl 


J 


20 


1.5" 


2 


CAl 


H 


S2 


1.5"-2.5" 


20 


CAl 


K 


25 


1.5" 


2 


CA2 


H 


22 


2"-3" 


18 


CAS 


J 


5 


1.8" 


4 


CAS 


H 


5 


1.8" 


4 


CAS 


K 


6 


1.8" 


4 


CA4 


H 


SO 


2.8" 


4 


El 


K' 


20 


1.5" 


5 


E2 


K' 


12 


1.5" 


4 



Notes: 
Col. (1) 
Col. (2) 
Col. (S) 
Col. (5) 



Observing run; abbreviation see text. 
Filter. 

Typical integration time on source in minutes. 
Number of observed galaxies; some galaxies were ob- 



served twice. 

Data reduction is performed with standard techniques 
for NIR images: dark subtraction, sky background sub- 
traction, bad pixel correction, and flatfielding (CAl-4: 
sky fiats, El-l-2: dome flats) (see e.g. MAGIC Observer's 
Guide [] or ESO-MIDAS User Guide, Volume B g). The 
sky background frames are determined by offsets, after 
each object frame, integrating sky images with the same 
time as used for the object images. Combining the indi- 
vidually reduced object frames leads to the flnal image 
of the galaxy. For the photometric calibration of the im- 
ages infrared standard stars (Elias et al. 
al. 



1982 



Persson et 

19981) are used giving a photometric error of A ~ 0.3 



mag. The data obtained at TIRGO are calibrated using 
published integrated aperture data. 

3. Analysis of the data 



Following Wakamatsu & Hamabe (1984) "bumps", which 
are symmetric to the minor axis, in cuts along or parallel 
to the major axis of edge-on disk galaxies can be asso- 
ciated with bars. This is indicated by the fact that the 



^ http:/ /www. caha.es/CAHA/Instruments 
^ http:/ /www. eso.org/projects/esomidas 
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bump becomes more and more insignificant at larger dis- 
tances from the plane (D'Onofrio et al. 1999| ). Using this 
technique bars in edge-on galaxies NGC 4762 (Wakamatsu 
fc Hamabe |l98l, NGC 1381 (de Carvalho & da Costa 



1983), ^CIC 5170 (Dett mar fc Barteldrees |l990b|), NGC 
4302 (Dettmar & Ferrara p^ ) NG C 325GB, NGC 5047, 
NGC 6771, IC 4767 (Liitticke |1999D, a nd at a marginal 
level in NGC 128 (D'Onofrio et al. |1999|) have been found. 



However, for edge-on galaxies these radial cuts are obvi- 
ously hampered by dust near the plane. Therefore previous 
studies using optical observations are limited to early type 
galaxies with low dust content or the resulting evidence 
for a bar remain uncertain (NGC 5170). The influence 
of dust is significantly reduced in the NIR enabling the 
identific ation of bars also in late type spirals (Dettmar & 
Ferrara [l99^ ). 

Therefore we have searched for bumps, which are sig- 
natures of a bar, in radial cuts of the galaxies in our NIR- 
sample. Additionally, the shape of the bumps in the cuts 
along and parallel to the major axis can be divided in 
profiles with a flat part and in profiles decreasing from 
the center of the galaxy throughout. For these purposes 
we have determined the surface brightness at the radial 
position of the top of the bump (i?tb) and at the end of 
the central bulge (i?cBu)- We call the feature in the cuts 
"strong bar signature" , if the gradient 

M(^tb) - Ai(-RcBu) 



bump 



-Rtb — ^CBU 



is smaller than 0.015 mag arcsec" -1 (Fig. 0and|). Other- 
wise the bump is called "weak bar signature" (Fig. ||). The 
error of this gradient is ±0.003 mag arcsec"^. An example 
for a galaxy without a bar signature is given in Figure ||. 

For galaxies showing a bar signature several charac- 
teristic quantitative parameters are derived from the cuts 
parallel to the major axis (for illustration see Fig. |l]): 

— BAL: The projected bar length is marked in the cut 
along the major axis by an increasing light distribution 
towards the center compared to the radial exponential 



light distribution (Pohlen et al. 2000) of the disk. 

— BAT: The bar thickness is defined by the vertical ex- 
tent of cuts showing the bar signature. 

— CBU: We introduced the length of the central bulge 
as the length of the structure being brighter than the 
bar signature. 

— BUL: The bulge length is marked by the increasing 
light distribution over the exponential disk well above 
the bar. 

— BPL: The length of the b/p structure is distance be- 
tween the maxima of the b/p distortion. This param- 
eter is obviously only deflnable for peanut bulges. 

However, taking into account seeing effects the measured 
values are only upper limits. 

For 17 galaxies of our sample it is not possible to search 
for bar signatures. Sometimes the heavy dust extinction 



NGC 2654 




-20 20 

R [arcsec] 



Fig. 1. Top: The peanut bulge of NGC 2654 in the NIR. 
CAl, 25min in H . Bottom: In cuts along and parallel to 
the major axis of NGC 2654 the parameter BAL, BAT, 
CBU, BUL, and BPL can be determined as explained in 
the text. The bar signature can be detected in cuts of 
the range z [—3", -1-3"]. The flat bump has a gradient 
Vbump = (16.90 - 16.75) mag *(34 - 18)-iarcsec-i = 0.01 
mag arcsec^ ^ (strong bar signature). The feature at the 
top of the bump at R = +34" can be associated with a 
spiral arm. 



near the plane, even in the NIR, prevents a final analysis. 
In other galaxies extra features (e.g. stars or possible spi- 
ral arms), asymmetrical structures, or low signal-to-noise 
ratios would lead to a too uncertain verification of a bar. 
Therefore the resulting sample is reduced to 43 galaxies. 

4. Results 

The derived bar signatures and parameters for our sample 
are presented in Table ||. Furthermore, we have computed 
characteristic ratios of the parameters (Table H) . 

In Paper I we have visually classified bulge shapes 
in three types of b/p bulges (1: peanut-shaped, 2: box- 
shaped, 3: nearly box-shaped) and one elliptical type (4). 
Regarding the ratio of the projected bar length to the 
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Fig. 2. Top: Contour plot of NGC 7332. CA4, 20min in H. 
Bottom: Cuts along and parallel to the major axis reveal a 
decreasing bump: Vbump ^ 0.03 mag arcsec"^ (weak bar 
signature). 



bulge length (BAL/BUL) we measure a clear trend with 
the bulge type (Fig. ^: The more prominent the b/p type, 
the larger the ratio. The mean value for type 1 bulges is 
BAL/BUL = 1.9 ±0.3, type 2: 1.6 ±0.2, type 3: 1.4 ±0.2, 
and type 4: 1.1 ± 0.1. However, we are dealing with small 
number statistics, N = 21, and the scatter is relatively 
high. This could be explained by the different morpholog- 
ical types contributing to these statistics, because early 
type galaxies have smaller ratios tha n late types (Fig. 
Additionall y, acc ording to Martin ( 1995 ) , Elm egreen & 
Elmegreen (1995), and Elmegreen et al. (1996) the rela- 
tive bar length {Rhar / R25) has also a high scatter pointing 
to a large range of bar lengths. 

Two other characteristic ratios, CBU/BPL and 
BAL/BPL, are only related to peanut bulges. They can 
be compared with predictions of N-body simulations and 
thereby test the evolution scenario of b/p bulges proposed 
in these simulations (see next section). The mean value of 
BAL/BPL is 2.7 ± 0.3 (Table |). A surprising result for 
CBU/BPL is the mean value of 1.0 with the small scatter 
of 0.1 (Tableland Fig. p. 




Fig. 3. Top: Contour plot of UGC 3354. CAl, 35min in 
H. Bottom: There is no bar signature visible in the cuts 
along and parallel to the major axis. 



Five out of the six galaxies without a bar signature 
(Table ||), but with a b/p bulge are thick boxy bulges 
(TBBs) (Dettmar & Liitticke |1999D. This group of galax- 
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Fig. 4. The bulge type (in the NIR) seems to be correlated 
with the ratio of the projected length of the bar (BAL) and 
the bulge length (BUL). Crosses: SO-Sab; squares: Sb-Scd. 
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Table 3. Connection between b/p bulges and bar signa- 
tures 



classified 


number 


of ga 


axies with 


bulge type 


no sig. 


WB 


' SB 2 


1 








8 


2 


0(1) 


2 


5 


3 


1 (4) 


5 


1 


4 


13 


3 






Notes: 

galaxies in brackets are excluded (see text) 
^: Weak bar signature 
^: Strong bar signature 



ies shows remarkable morphological features such as large 
asymmetries and substructures and the box-shape is more 
prominent in the outer part of the bulge, possibly pointing 
to a different formation scenario. Therefore this group of 
b/p bulges will be studied separately in a forthcoming pa- 
per (Paper III). The remaining galaxy (NGC3762) with- 
out bar signature and b/p bulge is classified as a weak 
type 3 bulge on an optical image (Paper I). Excluding 
the five TBBs from Table ^ there is an almost one-to- 
one correspondence between b/p bulges and bar signa- 
tures. The fraction of b/p bulges with such a signature 
is 95 % whereas over 80 % of the non-b/p bulges show no 
bar signature. Additionally, it should be mentioned that 
the three galaxies (NGC4026, NGC4762, and UGC9858) 
with a bar signature and without a b/p bulge have the 
smallest BAL/BUL ratios in our sample (Table ^ and Fig. 
^). Table ^ reveals a significant correlation on one side 
between strong bar signatures and prominent b/p bulges 
(type 1-1-2) and on the other side between weak bar sig- 
natures and less pronounced b/p bulges of type 3. This 
strong dependence is shown on a 1.0% significance-level 
by a x^-test. We do not find a correlation between the 
Hubble type and the weak or strong bar signature. 

For NGC 1175 (Fig. |) and NGC3203 we measure a 
strong bar signature. Therefore the classification as un- 
barred systems within the RC3 can be replaced. Since the 
bar classification for edge-on galaxies is difficult and lead 
to misclassifications in galaxy catalogues based on photo- 
graphic plates, our method of bar detection points at a 
solution for this problem. 

5. N-body simulation 

5.1. Model 

The origin of b/p bulges based on the evolution of bars is 
shown in several N-body simulations (e.g. Combes et al. 
199C, Pfenniger & Friedli 1991). To check the consistency 

we 



of simulations with the results of our observations, 
have investigated a representative simulation of a barred 
galaxy. For this purpose a particle position file of the final 
stage of a N-body simulation of a barred galaxy in three di- 




20 
R [arcsec] 



Fig. 5. Cuts along and parallel to the major axis of NGC 
1175 show strong bar signatures. This points to the exis- 
tence of a bar, which is not listed in the RC3. CA3, 6 min 
in H. 



mensions was kindly made available by Daniel Pfenniger. 
The simulation uses a particle-mesh code with 200.000 
particles described in Pfenniger & Friedli ( 1991 ). They 
choose axisymmetric initial conditions which are likely to 
develop a bar instability on their own and an initial mass 
distribution corresponding to a superposition of two ax- 
isymmetric Miyamoto-Nagai disks. The provided binary 
file represents the final time step of their simulation at T = 
5000 (~ 5 Gyrs) in which a bar is formed. The length unit 
is chosen to be one kpc (Pfenniger & Friedli 1991). For the 
transformation of the three-dimensional mass distribution 
of the particles to a two-dimensional surface brightness 
distribution we assume a constant M/L ratio and take 
into account the inclination (j) of the galaxy and the as- 
pect angle (0) of the bar. An edge-on bar is denoted with 
6 — 90°, while 6 — 0° refers to an end-on bar. 



5.2. Classification of the bulges 

In the investigated N-body simulation the peanut shape 
of bulges can be detected down to a galaxy inclination 
of ^^75°, if the bar is edge-on. If the galaxy is exactly 
edge-on, peanut bulges can be seen down to an aspect an- 
gle of the bar of ~ 77° . In general the simulation shows 
that the degree of boxiness of the bulge, reflected by the 
bulge type, is weakened with decreasing inclination of the 
galaxy (Fig. ||, Table At inclinations lower than ^ 75° 
it is not possible to distinguish between structures per- 
pendicular to the galactic plane and the bar itself. This 



value was also previously found by Shaw et al. (1990) as 
limit for the observation of b/p bulges. Regarding galaxy 
inclinations above this limit the b/p morphology is de- 
tected in the simulation whenever (j) is larger than ^^20°. 
In this way ~77% of projection angl es of bars result in 
b/p bulges (Table |). Combes et al. (|l990D find in their 
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simulations a similar percentage (85 %) for "favourable" 
projection angles to see b/p bulges in galaxies with an 
inclination between 80° and 90°. 

Since the frequencies of the simulated bulge types (Ta- 
ble ^ Col. 4) are derived only from a barred galaxy, these 
frequencies must be modified to compare them with our 
statistics of bulge types (Paper I) derived from all edge- 
on galaxies, with bar and without. Taking this into ac- 
count the percentages for the bulge types of the simulated 
galaxy have to be multiplied with the observed fraction 
of barred galaxies. Adding afterwards the fraction of un- 
barred galaxies to the percentage of bulge type 4 — N- 
body simulations (Pfenniger & Friedli 199lD suggest that 
unbarred galaxies have such a bulge type — leads to the 
modified frequencies of the bulge types (Table ^, Col. 5 
and 6). The errors of these frequencies depend on the un- 
certainties of the limits in the aspect angle of the bar {</>) 
(up to ±4°) and have values of ^2.5% for type 1 and 4 
bulges and '--^ 5 % for type 2 and 3 bulges. The observed fre- 
quencies of the bulge types (Paper I; type 1: 4.1 %, type 2: 
15.7 %, type 3: 25.2 %, and type 4: 55.0 %) are very similar 
(inside the errors) to the modified frequencies computed 
with a fraction of 55 % barred galaxies. This percentage is 
re vealed by studies fro m op tical surveys (e.g. RC3) (Ho et 
al. 1997 , Knapen et al. 200C, Paper I). However, the higher 
percenta ge of 70% derived fr om N IR studies (Mulchaey 
& Regan |l997l Knapen et al. |200C| , Eskridge et al. ^OOOj) 
seems to be more reliable, since the NIR-light traces the 
dominant, old stellar population and is relatively unaf- 
fected by dust. Using this fraction for computing the mod- 
ified frequencies, there is a difference of ^--^10% between 
observed and modified frequencies of galaxies with b/p 
bulges (type 1-3). There are less b/p bulges observed as 
computed from the investigated simulation coupled with 
NIR observations. 

The existence of galaxies with non-b/p bulges, but 
with bars (Tab. |^) are consistent with the simulation 
where almost end-on orientated bars lead to elliptical 
bulges. 

5.3. Quantitative parameters 

Along the bar major axis of the simulated galaxy, mea- 
sured in the face-on case (i = 0°), the decrease of the sur- 
face density is exponential (Pfenniger & Friedli 1991) and 
the length of the bar can be determined to ^ 20 kpc. This 
length exactly fits inside the radius of cor otation which 
occurs at 10 kpc (Pfenniger & Friedli 1991). The axis ra- 
tio of the bar is b/a = 0.7 ± 0.1 (Fig. y, top right) and is 
inside the range of observed values (0.2 - 1) (Martin [l995| ) . 

Regarding the galaxy edge-on {i = 90°) with an edge- 
on orientated bar {(j> = 90°) is the best way to derive the 
quantitave bar parameters (introduced in section 3) of the 
simulated galaxy, because bar and b/p structures are most 
prominent at this geometry. However, the cuts along and 
parallel to the major axis do not reveal the typical bar 



Table 4. Bulge types in the N-body simulation 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


bulge 


i 


<^ 


fre- 


mod. 55 % 


mod. 70% 


type 






quency 


frequency 


frequency 


1 


90° 


77° -90° 


14.4 % 


7.9 % 


10.1 % 


2 


90° 


47°-76° 


33.3 % 


18.3 % 


23.3 % 


3 


90° 


16°-46° 


34.5 % 


19.0 % 


24.1 % 


4 


90° 


0°-15° 


17.8 % 


54.8 % 


42.5 % 


1 


85° 


81°-90° 


10.0 % 


5.5 % 


7.0 % 


2 


85° 


52°-80° 


32.2 % 


17.7 % 


22.5 % 


3 


85° 


18°-51° 


37.8 % 


20.8 % 


26.5 % 


4 


85° 


0°-17° 


20.0 % 


56.0 % 


44.0 % 


1 


80° 


85°-90° 


5.6 % 


3.1 % 


3.9 % 


2 


80° 


57°-84° 


31.1 % 


17.1 % 


21.8 % 


3 


80° 


20°-56° 


41.1 % 


22.6 % 


28.8 % 


4 


80° 


0°-19° 


22.2 % 


57.2 % 


45.5 % 


1 


75° 


89°-90° 


1.1 % 


0.6 % 


0.8 % 


2 


75° 


66°-88° 


25.6 % 


14.1 % 


17.9 % 


3 


75° 


30°-65° 


40.0 % 


22.0 % 


28.0 % 


4 


75° 


0°-29° 


33.3 % 


63.3 % 


53.3 % 



Notes: 

Col. (1): Bulge type as defined in Paper I. 

Col. (2): Inclination of the simulated galaxy. 

Col. (3): Aspect angle of the bar (90° — edge-on and 0° — 

end-on). 

Col. (4): Frequency of the bulge type of the simulated galaxy, 
viewed under different aspect angles, at a constant inclination. 
Col. (5) : Modified frequency of the bulge type of the simulated 
galaxy, viewed under different aspect angles, at a constant 
inclination taking into account that only 55 % of all galaxies 
are barred and the remaining galaxies without a bar (45 %) 
have a bulge of type 4. 

Col. (6): As Col. (5), but with 70% barred and 30% unbarred 
galaxies 



signature (Fig. ^ which is visible in the NIR observa- 
tions. The three components bulge, bar, and central part 
of the bulge cannot be distinguished as it is possible in 
the observations (cf. Fig. ^ ^ and ||). Due to the lack of 
a bar signature the bar thickness as well as the central 
bulge length (CBU) cannot be determined. The ratio of 
BAL/BUL = 1.25 ± 0.2, taking into account BAL = 20 
kpc and BUL = 16 kpc (Fig. is not consistent with 
the observed ratios for peanut bulges (1.9 ± 0.3) (Fig. 
while BAL/BPL = 2.5±0.2 (BPL -8 kpc; Fig. |) is con- 
sistent with observed ratios in peanut bulges (BAL/BPL 
= 2.7 ± 0.3). While in observations the positions of the 
maxima of the b/p distortion coincide with the limits of 
the central bulge, in the N-body simulation the positions 
of the maxima are not related to any other features in the 
profiles. 



8 



R. Liitticke et al.: Box- and peanut-shaped bulges. II 




published simulations of barred galaxies (e.g. Pfenniger 



1999) to get more relevant comparisons with observations. 



Fig. 7. Cuts along and parallel to the major axis of the 
simulated galaxy with i — 90° and cj) = 90°. Cuts with 
the same distance from the major axis are averaged. BPL 
= 8 kpc and BUL — 16 kpc can be derived from the cuts, 
while a bar signature is not detectable. BAL = 20 kpc is 
determined at i — 0°, see Fig. 2 in Pfenniger & Friedli 



(1991) 



5.4- Relevance of the simulated galaxy for observations 

The differences between the galaxy derived from the N- 
body simulation and the observed galaxies likely result 
from the initial conditions of the simulated model. Pfen- 



niger & Friedli ( |l99l[) are just using two merged disks 
with the same rotation properties instead of an indepen- 
dent bulge component in their simulation. The inner disk 
represents the "bulge" and is chosen with the main pur- 
pose "to have sufficient freedom to shape the initial rota- 
tion curve" . Therefore the initial resemblance with a real 
bulge is only approximate. Using such a bulge component 
in the simulation could also be the reason for the lack of 
the bar signature when the galaxy is seen edge-on, since 
the inner disk could cover the signature. Another expla- 
nations for the fact that no bar signature is visible could 
lay in the limited spatial resolution or in the lack of gas 
in the simulation, which influences the bar structure. It 
is also possible that the exponential decrease of the bar 
surface brightness and the smooth end of the bar prevent 
the detection of a bar signature. 

Additionally, the investigated simulated galaxy is only 
a single, although representative, example of a barred 
galaxy. Variations of the initial conditions such as the 
bulge-to-disk ratio, Q parameter, gas fraction and dis- 
tribution, and halo structure end up in different bar 
strenghts, shapes, thickness, lenghts, and density proflles 



in N-body simulations (e.g. Friedli & Benz 1993). There- 
fore conclusions drawn from comparisons between our 
statistics of observed bulge types (Paper I) and results 
from the simulated galaxy, which is investigated in this 
study, are limited. Future studies have to analyse recently 



6. Discussion 

6.1. Bars and b/p bulges 

The strong correlation of the bar signature with b/p 
bulges in our study shows that b/p bulges are well ex- 
plained by dynamical processes in bar potentials. The ra- 
tio of the bar length to the length of the b/p distortion 
(BPL) (2.5 ± 0.2) in the simulated barred galaxy and the 
observed ratio of BAL/BPL (2.7 ± 0.3) in galaxies with 
peanut bulges is consistent with such a scenario. 

The correlation of the bar signature with b/p bulges 
is previously claimed, using the same detection method, 
by the investigations of individual galaxies, namely NGC 
1381 (de Carv alho fc da Costa ^98^, NGC 5170 (Dettmar 



Ba rteldrees |l990bD , and NGC 4302 (Dettmar & Ferrara 
1996|) . Many further evidences for the connection between 



bars and b/ p bulg es are reported in th e liter ature: Bettoni 
& Galletta ( |l994 and Quillen et al. ( |l997| ) find them in 
their investigation of the intermediately inclined galaxies 
NGC 4442 resp. NGC 75 82. F urthermore, N-bo dy simula- 
tions (e.g. Combes ct al. 1990, Raha et al. 1991 , Pfenniger 
& Friedli 1991 ) show that the evolution of a b/p bulge is 
coupled with the development of a bar. 

Independent evidence is obtained by observing the gas 
kinematics in Nil and Ha. The characteristic "figure-of- 
eight" rotation curve which is a strong signature for the 
presence of a bar and is detected in many b/p bulges: e.g. 
NGC 5746, NGC 5965 (Kuijken & Merrifield |l99|) , NGC 
2683 (Merrifield 1996), and additional ten further galax- 
ies by Bureau & Freeman ( 1999D and five by Merrifield & 
Kuijken ( |1999 ). Control samples of non-b/p bulges show 
no k inematical signature of a bar (Bureau & Freeman 



1999| , Merrifield & Kuijken [19991) . In two b/p bulges (ESQ 



597- 36 and NGC 33 90) w ithout kinematic bar signature 
Bureau & Freeman ( 199£ ) find indications for accretion 
events. However, since the positioning of the slit was diffi- 
cult (Bureau & Freeman 1999), dust in the galactic plane 
could affected their observations. Therefore bars can not 
be excluded in these two galaxies. A merger- induced bar 
formation with a resulting b/p bulge (Mihos et al. 1995) 
could be a possible scenario for these galaxies. The four de- 
tected non-b/p bulges: NGC 5907 (Miller fc R ubin |1995D , 
IC 5096, ESQ 240- 11 (Bure au fc Freeman |1999| ) , and NGC 
681 (Merrifield & Kuijken 1999) with a bar signature in 
the gas kinematics can be explained by our analysis of the 
N-body simulation with bars which arc nearly end-on (see 
next subsection). 

Signatures of the gas kinematics in barred potentials 
are also detected in retrograde orbits of NGC 128 (Em- 



sellem & Arnault 1997), the galaxy wi th on e of the most 



prominent b/p bulges (D'Onofrio et al. 199£ ). Direct kine- 
matic evidence for bar streaming motions in two galaxies 
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with b/p bulges ( NGC 3079 and NGC 4388) is presented 
by Veilleux et al. ( |1999D . Furthermore, in HI the presence 
of a bar is also verified in one galaxy with a peanut bulge 
analysing the position velocity diagram (IC 2531: Bureau 
& Freeman 1997). The gas kinematics in the central region 
obtained in C O (J = 2-1 line) in the case of NGC 3628 
(Ren ter et al. p^9l| ), NGC 891 (Garcia-Buri Uo fc Guehn 
1995| ), and NGC 4013 (Garcia-BuriUo et al. |1999D again 



point to the presence of a bar. 

The conclusion (often given in the literature, e.g. Pfen- 
niger & Friedli 1991), that b/p bulges are thick bars, can- 
not be supported. From the observational point of view 
bars and bulges, if b/p-shaped or not, are obviously in- 
dividual components with different lengths (e.g. Fig. |l|). 
Stellar orbits in a barred potential build up b/p structures 
and the bar itself. Therefore b/p bulges are the composi- 
tion of a spheriodal bulge and a b/p structure representing 
the thick part of the bar. 

The non-detection of bars (Table ^ in all observed 
TBBs (~ 5 % of all b/p bulges belong to this group. Paper 
III) is an indication for the hypothesis that the boxiness 
of these bulges is not related to the the presence of a bar. 
This idea gets more convincing support from the fact that 
the extent of the box shape is too large to result from 
an usual bar potential (Paper III). However, the lack of 
a bar signature does not necessarily mean that there is 
no bar, since weak bars, which have exponential surface 
brightness profiles (Elmegreen & Elmegreen 1989), could 



remain undetecable by the here presented search for a bar 
signature. 

In this way also the missing bar signature in the sim- 
ulated barred galaxy is explainable, because the surface 
density distribution along the simulated bar is exponential 
and additionally, the bar axis ratio of this galaxy points 
to a weak bar {b/a > 0.6, Martinet & Friedh |1997D. 



6. 2. Explanations for various bulge shapes 

The connection between b/p bulges and bars can be un- 
derstood as interplay between resonant bending fed by 
vertical diffusion of orbits and instabilities similar to the 



fire-hose one (Bureau 199S, Friedli 1999). While there is 
evidence for this connection in the most galaxies with b/p 
bulges, the question for the different shape of the b/p 
bulges remains to be answered. 

A correlation between a strong signature of a bar, 
marked by constant surface brightness in radial direction 
(cf. section 3), and prominent b/p bulges is found in our 
study. Thereby the strong bar signature can be associated 
with a flat bar, which is known to be strong (Elmegreen 
& Elmegreen 19891 ). Therefore the correlation can be in- 
terpreted by the assumption that strong barred galaxies 
are responsible for prominent b/p bulges, while galaxies 
with a weak bar have bulges with a lower degree of box- 
iness. N-body simulations by Friedli (1999) support this 
explanation for different bulge shapes. In this way it is 



also possible to explain partly the observed correlation of 
BAL/BUL with the bulge type (cf. section 4) by larger 
bars in galaxies with prominent b/p bulges, because flat 



bars tend to be larger (Elmegreen & Elmegreen 1995, 
Elmegreen et al. 1996). Since this correlation holds for 
early as well as for late types, it cannot be explained by 
the fact that the same bar will produce more pronounced 
b/p bulges, the smaller the pre-existing spheroidal bulge 
is. However, having this in mind it is not surprising that 
the weak bar of the here investigated simulated galaxy can 
produce a peanut bulge. Furthermore, one could expect 
that peanut bulges more often occur in late type spirals, 
but the observed fraction of peanut bulges in galaxies in 
the range of SO-Scd is rather constant and equal to zero in 
spirals later than Sc (Paper I). This result could be inter- 
preted by a dependence of the bar strength on the Hubble 
type. Such a correlation of bars is found by Elmegreen & 
Elmegreen ( |1985D and Elmegreen e t al. ( |l996[ ). However, 
in contrast to them Seigar & James ( |1998| ) do not find that 
flat (strong) bars are more common in early Hubble types 
and exponential (weak) bars in later types. This state- 
ment is in agreement with our results for the profiles of 
bar signature concerning the Hubble type, but it is uncer- 
tain, if their results concerning bar morphology detected 
in face-on galaxies can be transformed to bar signatures 
in edge-on galaxies. 

Additionally, for the explanation of bulge shapes it has 
to be taken into account that gas mass and a central black 
hole effect the shape of b ulges . A significant fraction of gas 
mass (Berentzen et al. 1998| ) or a supermassiv e cent ral 
black hole can weaken the b/p structure (Friedli 1999|) . 

Explaining bulge shapes by aspect angles of bars was 
introduced by N-body simulations (Combes et al 



1990 



Pfenniger & Friedli 1991). Kinematical bar diagnostics in 
edge-on spiral galaxies using simulations of families of pe- 
riodic orbits and hydrodynamical simulations are also con- 
sistent with observations. They confirm the dependence 
of the boxiness on the aspect angle of the bar (Bureau 
Athanassoula & Bureau 1999). An 



Athanassoula 1999 



exact determination of aspect angles is not possible, be- 
cause the envelope of the signature of the outer bar region 
in the posi tion- v elocity diagram is too faint (Athanassoula 
& Bureau |l999| ). 

The strong correlation of BAL/BUL with the bulge 
type and the agreement of the observed ratio BAL/BUL 
only of the most prominent b/p bulges, 1.9 ± 0.3, with 
that ratio obtained for face-on galaxies, 2.3±0.7 (Athanas- 
sou la fc Martinet |1980| , Baumgart & Peterson |1986| , Mar- 



tin 11995), shows that the differences in BAL/BUL can be 
interpreted by different aspect angles of bars. Therefore 
our analysis points also to a dependence of the shape of 
a bulge on the aspect angle of the bar. In this way and 
consistent with our investigation of the simulated barred 
galaxy peanut bulges result from bars seen nearly edge-on 
(large projected bar length, BAL), boxy bulges from in- 
termediately inclined bars, and elliptical ones from end-on 
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bars (small BAL) or they have no bar at all. Thereby bars 
detected in non-b/p bulges can be explained as well as 
differences between the fraction of barred galaxies (70 % 
in the NIR) and galaxies with b/p bulges (45 %), although 
also other parameters might play a role (weak bar, large 
gas mass, or dominant bulge component). The dependence 
of the bulge shape on the aspect angle of the bar is also 
supported in our study by the modified frequencies of the 
bulge types derived in the simulation. These frequencies 
are consistent with the observed frequencies of the bulge 
types, if one keeps in mind that only a single simulated 
barred galaxy is compared with all morphological types 
of disk galaxies. Larger bulge components in the simula- 
tion would weaken the b/p structures and thereby reduce 
the modified frequency of b/p bulges. Stronger bars would 
have the contrary effect on the modified frequency. Addi- 
tionally, many galaxies of our investigated sample with 
peanut bulge and strong bar signature show features at 
the ends of the bumps, which can be associated with spiral 
arms extending outward from the ends of a bar (Baggett 
& MacKenty 1996 ). Therefore it can be suggested that 
bars of these galaxies are indeed nearly edge-on. 

There exists only a few rough estimates of aspect an- 
gles of bars in edg e-on g alaxies from observations. Waka- 
matsu & Hamabe { \l984 ) adopt for NGC 4762 (bulge type 
3 in the optical, type 4 in the NIR; Paper I) an aspect 
angle of the bar oi (p — 33° ± 10° (although a lower 
limit of 21° is also possibl e from their argumentation) and 
Garcia-BuriUo & Guelin ( |1995D for NGC 891 (bulge type 
2) ^ 45°. These results are in agreement with a depen- 
dence of the bulge shape on the bar aspect angle. 

6.3. A new correlation 

The coincidence of the positions of the maxima of the b/p 
distortion with the limits of the central bulge, which is 
found here for the first time, cannot finally be explained 
(CBU/BPL= 1.0 ± 0.1). Selection effects due to different 
aspect angles of the bar, which would infiuence the length 
of the b/p distortion (BPL), can be neglected, since we 
have shown that peanut bulges, for which BPL is defined 
only, can be associated with nearly edge-on bars. The co- 
incidence points also to a possible connection between the 
bar length and the length of the central bulge (CBU), since 
BPL is connected with BAL. 

Therefore possible explanations could be related to sec- 
ondary bars (e.g. FriedU |1999| , Greusard et al. |000|). This 
idea is supported by the result from N-body simulations 
(Friedli & Martinet 1993| ) that the corotation of the inner 
bar roughly coincides with the location of the horizontal 
inner Lindblad resonance (ILR). At this position also the 
maximum of the elevation of the peanut structure occurs, 
because the vertical ILR is very close to the horizontal 
ILR and the conjunction of these two resona nces build 
up the peanut structure (Combes et al. 1990). However, 
the observed ratio of the bar lengths of the primary to 



the secondary, which ends near corotation (Greusard et 
al. |200C| ), varies between 4.0 and 13.4 (Fr iedli e t al. |1996D , 
respectively 2.7 and 7.1 (Greusard et al. 2000| ), while the 
ratio of bar length to length of the central bulge is much 
smaller and varies only between 2.4 and 3.1. 

The central bulge could also be associated with a pri- 
mordial bulge, which influences the bar structure and the 
position of the horizontal ILR. Since at this location also 
the peanut structures is most prominent, the extent of 
the b/p distortion would depend on the mass and length 
of the primordial bulge. Such a hypothesis is supported by 
N-body simulations in which the size of the evolving bar is 
direc tly related to the initial size of the bulge (Pfcnniger 
1999 ). In this way the nearly constant ratio of bar length 
or BPL to CBU would be a natural result. 

The correlation between CBU and BPL and especially 
the CBU/BPL ratio of 1 could also be explained by nu- 
clear rings. They corres pond to the horiz ontal ILR of bars 
(Buta & Crocker |1993| , Buta & Combes |1996| ) and there- 
fore occur at the same radial position as the maximum of 
the elevation of the peanut structure. A nuclear ring could 
be responsible for the increase of the surface brightness, 
which we observe in radial cuts near the plane (Fig. ^ 
H, 1^), or the ring is only the border for the central bulge. 
However, we have not seen any structure inside the central 
bulge pointing to rings, but it is likely that our resolution 
is too low to detect them. 

6.4- Bar thickness 

Since we know now with high probability the fact that 
bars in peanut bulges are edge-on orientated, it is for the 
first time possible to determine the thickness of bars by 
observations. From the mean value of BAL/BAT for the 
galaxies with peanut bulges a ratio of bar length to thick- 
ness of 14 ± 4 (taking the error in BAL and BAT into 
account) is derived. This value is at the upper limit of the 
values estimated in the literature (ranging from 4 to 15, 
Sellwood & Wilkinson 1992). It must be mentioned that 
bars are more extended in the NIR than in the optical 



(Friedli et al. 1996, Chapelon et al. 1999). However, from 



a sample of nine barred galaxies (Wozniak et al. 1995 



Friedli et al. 199(;) investigated in the /- and K-hand we 
have computed that this effect is only of the order of ~ 5%. 
Therefore our rather large value for the ratio is not an ef- 
fect of the observed wavelengths. 

7. Conclusion 

We have analysed a sample of edge-on disk galaxies in the 
NIR and images derived from a N-body simulation of a 
barred galaxy using cuts parallel to the major axis and 
our classification of bulge shapes introduced in Paper I. 
Our main conclusions are: 

— Cuts along and parallel to the major axis of edge-on 
galaxies show a strong correlation between bars and 
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b/p bulges. This correlation is for the first time verified 
in a larger sample of galaxies (N = 60). 

— The investigation of the bar signatures in the NIR re- 
veals a correlation between prominent b/p bulges and 
strong bar signatures pointing to a dependence of the 
boxiness of a galaxy on the bar strength. 

— The ratio of the projected bar length (BAL) to the 
length of the bulge (BUL) correlates with the bulge 
types, i.e. with the bulge shape. This correlation is the 
first evidence from surface photometry for an explana- 
tion of the different degree of boxiness by the aspect 
angle of a bar. The maximum of this ratio for galaxies 
with peanut bulges is consistent with ratios observed 
in face-on galaxies. 

— The ratio of BAL to the length of the b/p distortion 
(BPL) has a mean value for the observed peanut bulges 
of 2.7 ± 0.3. This value is consistent with the investi- 
gated N-body simulation where a ratio of 2.5 ±0.2 can 
be derived for an edge-on galaxy with an edge-on bar. 

— The frequencies of bulge types in the simulated galaxy, 
seen under different aspect angles of a bar and un- 
der high inclinations, are consistent with the observed 
frequencies corrected for the total fraction of barred 
galaxies. These results support the interpretation that 
bars viewed under different angles in edge-on galaxies 
play a major role for the bulge shape. However, the 
viewing angle is not the only parameter, since the size 
of the pre-existing spheroidal bulge and various bar 
shapes and types also influence the bulge shape. 

— Interpreting galaxies with peanut bulges as galaxies 
with a bar seen edge-on the ratio of bar length to thick- 
ness, BAL/BAT = 14 ± 4, can be directly measured for 
the first time. 

— From the observational point of view bars and b/p 
bulges are different components revealed by their dif- 
ferent lengths. The inner part of barred disk galaxies 
can be described by three components: the spheroidal 
bulge (including the central bulge), the bar (=thin 
bar), and the b/p structure (=thick bar). 

— We have found a new size relation of the positions of 
the maxima of the b/p distortion with the extent of the 
central bulge (CBU/BPL=1 ± 0.1). This size relation 
is most likely related to the in-plane ILR of the bar. 
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Table 2. Results of the structure analysis 
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Table 2 continued. 
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14 
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Notes: 

Col.(l): Galaxies are ordered by RA for each observing run separated by lines (TI, CAl, CA2, CA3, CA4, El, and E2 ). If a 
galaxy is observed twice, only the data with higher quality are listed. 

Col. (2): Bulge type as defined in Paper I in the NIR, except *: optical bulge type (Paper I), because signal-to- noise ratio in the 
NIR image is too low. 

Col. (3): Bar signatures: - = no signature, WB = weak bar signature, SB = strong bar signature, ? = uncertain (see text) 
Col. (4) - (7): Abbreviations see text. ? = BAL is larger than the size of the image. 
Errors of the measurements : BAL and BUL ~20%, BAT ~30%, CBU and BPL ~ 10%. 



